Renormalization of the BCS-BEC crossover by order parameter fluctuations 
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We use the functional renormalization group approach with partial bosonization in the particle- 
particle channel to study the effect of order parameter fluctuations on the BCS-BEC crossover of 
superfluid fermions in three dimensions. Our approach is based on a new truncation of the vertex 
expansion where the renormalization group flow of bosonic two-point functions is closed by means 
of Dyson-Schwinger equations and the superfluid order parameter is related to the single particle 
gap via a Ward identity. We explicitly calculate the chemical potential, the single-particle gap, and 
the superfluid order parameter at the unitary point and compare our results with experiments and 
previous calculations. 
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The BCS-BEC crossover in a two-component Fermi 
gas has attracted the attention of theorists for several 
decades 

BiSSi- It is generally accepted that the na- 
ture of the superfluid state exhibits a smooth crossover as 
a function of the dimensionless parameter l/kpa Sl where 
kp is the Fermi momentum and a s is the s-wave scatter- 
ing length in vacuum. While for a small negative scat- 
tering length, i.e., l/kpa s <C — 1, the paired state gener- 
ated by the attractive interaction is a collection of spa- 
tially extended Cooper pairs (BCS limit), in the opposite 
limit l/kpa s 3> 1 the superfluid state can be viewed as a 
Bose-Einstein condensate of tightly bound fermion pairs 
(BEC limit). Of particular interest is the unitary point 
l/kpa s — where the scattering length diverges and kp 
sets the only length scale of the system. In this regime 
quantitative calculations are difficult because there is no 
small parameter to justify approximations. 

In the past few years several observables such as the 
chemical potential and the quasi-particle gap have been 
determined experimentally at the unitary point 
lioj ]. but there is still some uncertainty in the precise nu- 
merical values of these quantities. The unitary point has 
also been studied theoretically using Monte Carlo sim- 
ulations [ill, |l2j and various analytical methods based 
on field theoretical techniques [l^.Tl3. [TEL 16 L 17 1 or the 
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functional renormalization group (FRG) 1 



but also in this case the theoretical results have not con- 
verged yet. In such a situation it is desirable to study 
this problem using new approximation strategies which 
are complementary to previous calculations. In this work 
we shall therefore develop a novel FRG approach for the 
BCS-BEC crossover which is based on a suitable trunca- 
tion of the vertex expansion using skeleton equations and 
Ward identities. For fixed density we explicitly calculate 
the chemical potential, the single-particle gap, and the 
superfluid order-parameter at the unitary point in three 
dimensions, and compare our results with experiments 



and with previous calculations. 

We consider a system of neutral fermions with energy 
dispersion e/~ = fc 2 /(2m) and a short-range attractive 
two-body interaction g p depending on the total momen- 
tum p of a fermion pair. After decoupling the interaction 
in the particle-particle channel using a complex Hubbard- 
Stratonovich field x, the Euclidean action of our model 
can be written as S = So + Si, with Gaussian part So 
and interaction Si given by [22j] 



So 



Si = 



{-iu + £k)^Kvi>K<r + / 9 p 1 XpXp, (1) 



K 



[ipK+P^-KlXP + Ip-KllpK+PtXp]- ( 2 ) 



P JK 



Here, the energy = — \i is measured relative to 
the chemical potential [i, and the anti-commuting fields 
ipKa and ipKtr represent fermions with energy-momentum 
K = (ilu, fc) and spin projection a. The complex bosonic 
field xp is conjugate to the fluctuation of the superfluid 
order parameter with energy- momentum P = (iu>,p). 
For convenience we choose our sign convention such that 
g p > for attractive interactions. 

To derive FRG flow equations for our model we intro- 
duce a cutoff A into the Gaussian propagators appearing 
in Eq. ([T]) and consider the evolution of the generating 
functional of the irreducible vertices as the cutoff is re- 
duced [23|, [24| • The physical vertices are then recovered 
for A — > 0. The FRG equations for the irreducible ver- 
tices of the above mixed boson-fermion theory follow as 
a special case of the general FRG flow equations derived 
in Ref. [25(. In contrast to a previous FRG calculation 
[2H , we use here a scheme where the cutoff is introduced 
only in the bosonic part of the Gaussian propagator. The 
advantage of our boson cutoff scheme is that the initial 
condition for the fermionic self-energy is simply given 
by the self-consistent Hartree-Fock approximation (i.e. 
the BCS approximation), while the initial vertices in the 
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bosonic sector are given by closed fermion loops with an 
arbitrary number of bosonic external legs 25[ . In partic- 
ular, at the initial scale the bosonic two-point functions 
are given by the ladder approximation. In order to ob- 
tain numerically tractable FRG equations, we have to 
make further approximations: first of all, we neglect the 
momentum-frequency dependence of the vertices with 
one boson and two fermion legs, replacing these vertices 
by a momentum- and frequency-independent coupling 
7 = 7a- Moreover, we completely ignore vertices with 
two fermion legs and more than one boson leg. Within 
these approximations, the FRG flow equations for the 
anomalous (A) and normal (£) fermionic self-energies in 
our cutoff scheme are 



d A A(K) = 7 d A ( X ) 



7 



[Fp — Fp]A(P — K) , (3) 



which are shown graphically in Fig. [T] (a). Here, A(K) 
and B(K) are the anomalous and normal component of 
the fermionic single-particle Green function, which are 
related to the self-energies by 



A{K) = 



A(K) 



B{K) = 



D{K) ' 
G^j-K) 
D(K) 



(5) 
(6) 



where G -1 (iT) = iu - f fe - T,(K) and D(K) = 
G- 1 {K)G- 1 {-K) + \A(K)\ 2 . The functions F l j in 
Eqs. (J3]) and ([4]) are the bosonic single-scale propagators 
associated with longitudinal (upper index £) or transverse 
(upper index t) fluctuations of the field \- Choosing for 
simplicity a sharp momentum cutoff in the bosonc sector, 
the single-scale propagators are Fp = —S(A — \p\)Fp , 
where the bosonic propagators Fp can be expressed in 
terms of the irreducible bosonic self-energies IP-? (P) (po- 
larizations) as 



bp bp 
bp bp 



N(P) 



II" (P) 



n"(p) 



(7) 



with 



N(P) = L 9p ] 



^"(P)]^ 1 



n"(p)] + [n et (p)} 2 . (8) 



In order to determine the fermionic self-energies from 
Eqs. ([3]) and ([4| we need additional equations for the 
polarizations n* J '(P) and for the flowing order param- 
eter (x). Instead of explicitly writing down the FRG 
flow equations for II y (P), we shall use skeleton equa- 
tions (which follow from Dyson-Schwinger equations 
(2||) to relate the bosonic self-energies to the fermionic 
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n"(P) = 



kxp Vx-p ix-p 

\ 6 -26 + t> 

* K\J-K + P * K\J-K-P \Jk + ] 
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(c) 
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Z K\J-K + P * K\J-K-P 

KXP VX-P 



FIG. 1: (Color online) (a) Diagrammatic representation of 
our approximate FRG flow equations (J3j [4} for the fermionic 
self-energies A(K) and T,(K). (b) Exact skeleton equations 
for the bosonic self-energies U u (P) , E* (P) , It" (P) . (c) Exact 
FRG flow equation for the order parameter (x)- Solid arrows 
represent fermionic propagators, dashed arrows represent the 
complex boson fields \ an d X> dashed lines without arrow 
represent longitudinal components (p e of % — (<P £ + *¥ , *)/V / 2, 
while wavy lines represent the transverse components <^'. 



ones [27]. Graphically, the exact skeleton equations for 
the irreducible bosonic self-energies IP- 7 (P) are shown in 
Fig. [IJb) • Within our truncation where the three-legged 
boson-fermion vertex is approximated by a flowing cou- 
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pling 7 = 7a, these relations become 



n"(P) 



n"(P) 



Tl et (P) 



K 



K 



B{K)B{-K + P) - A( J PC)A(K + P) 



fP- 



-P) 



(9a) 



B{K)B{-K + P) + + P) 
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+ (P--P)J , 

B(K)B(-K + P) 



(9b) 



(9c) 



To close our system of flow equations, we still need an 
equation for the vertex renormalization factor 7 and the 
flowing order parameter {x) appearing in Eq. ([3|) . In our 
boson cutoff scheme, the flow of (x) is driven by the ir- 
reducible vertices with three bosonic legs, T m , T ut , r , 
as shown graphically in Fig. [ljc). The crucial point is 
now that from the requirement that the FRG flow pre- 
serves the gapless nature of the Bogoliubov-Anderson 
(BA) mode [26| we can easily obtain an expression for 
<9a (x) without explicitly considering the RG flow of the 
bosonic three-legged vertices. Therefore we note that the 
condition for the vanishing of the gap of the BA mode is 



9o 



-1 , n"(o) = 



(10) 



This condition is obviously satisfied at the initial RG 
scale A = Ao because the ladder approximation with 
self-consistent Hartree-Fock propagators is conserving. 
To make sure that the BA mode remains gapless for 
any value of A we simply require that Eq. (|10[) remains 
valid during the entire RG flow. With n 4t (0) given by 
Eq. (|9bj) . this is an implicit relation between 8aA(K) 
and d\Ti(K). By demanding that this relation is consis- 
tent with Eqs. and ((H) we can uniquely fix the RG 
flow of the order parameter (x)- Finally, using the U(l)- 
gauge symmetry of the action given in Eqs. ([1]) and ([2|) 
we can derive a Ward identity which relates the ratio of 
the anomalous self-energy and the superfluid order pa- 
rameter to the vertex renormalization factor 7, 



A(0) = 7 <X) 



(11) 



Eqs. (f3"))- (fTTj) form a closed system of integro-differential 
equations for the fermionic self-energies A(if), T,(K), the 
vertex renormalization factor 7 and the order parameter 
(x), which should be solved with the initial conditions 
A(K) Aa = (x)a = A and B(K) Ao = 0. Here, A is the 
single-particle gap in the BCS approximation. 

The numerical analysis of the system of coupled 
integro-differential equations ((3|)- (fTT|) is beyond the 
scope of this work. Here, we further simplify these equa- 
tions by neglecting the momentum-dependence of the 



fermionic self-energies and keeping only the linear fre- 
quency correction to the normal self-energy, replacing 
A(K) -> A and E(K) -> £ - (Z^ 1 - where 



Z = 



1 



1 - 



dS(V,0) 
d(iui) 



(12) 



j=0 



is the inverse flowing wave function renormalization fac- 
tor. The flowing single-particle propagators are then 
given by A(K) = ZA(K) and B(K) = ZB(K), where 



A(K) 
B(K) 



El) 



-A/(co 2 - 



El) 



(13) 
(14) 



with E k = [il + A 2 ] 1 / 2 , £ fc = e k - /}, l k = Ze k , p, = 
Z(fj,~H), and A = ZA. It should be noted that A can be 
identified with the physical single particle gap which can 
be measured, e.g., in tunnelling experiments. With these 
approximations, our system of flow equations reduces in 
D dimensions to 

AdAp. = W - 7 (A/k F ,o) D e_F,o 

g[F^ + P«-2iP«]P(P),(15) 



A9 A ln 7 = -(j/A)(A/k F , ) D e F , Q 

^-[Ff - Fp]A(P) . 



(16) 



The wave function renormalization factor is determined 
by Ad\Z = rjZ ' , with the flowing anomalous dimension 



>1 



7 (A/A: 



F,0) 



El 



2?<D£a 



2iFt 



El) 2 



(17) 



Due to the sharp momentum cutoff we may set P = 
(iu>, A). The FRG flow is further constrained by the con- 
dition (fit)]) that the BA mode is gapless and by the rela- 
tion 



A = Z 1 { x ) , 



(18) 



imposed by the Ward identity (|11|) . The dimension- 
less interaction terms Fp appearing above are defined 
by Fp = Z 2 ~fVoFp , and e^o, k F ,o, and vq denote the 
Fermi energy, the Fermi wave vector and the density of 
states at the Fermi energy (per spin projection) of a non- 
interacting system which has exactly the same density as 
our interacting system at the initial scale A = Ao. 

For convenience, we work with a momentum- 
independent bare coupling g p — > go- In dimensions 
D > 2 this gives rise to an ultraviolet divergence in the 
BCS gap equation which also appears in the polariza- 
tions ft"(P) = n l4 (P)/(Z 2 7^ ). For D > 2 we may ab- 
sorb this divergence into the bare coupling by introduc- 
ing the dressed coupling g via (Z 2 -fg)~ 1 — (Z 2 jgo)~ 1 — 



0, 







A/cf,o 



FIG. 2: (Color online) RG flow the three-legged vertex 7 and 
the wave function renormalization factor Z at the unitary 
point (kFa a = 00) in three dimensions. 



V^ 1 £ fc (2e fe ) _1 , so that the constraint (fTOj) that the BA 
mode remains gapless turns into 



1 



Z^g V 



1 

2E k 



1 

2^ 



(19) 



In D — 3 the dressed coupling is related to the s-wave 
scattering length a s via g — — 4ira s /m. It is intriguing 
to note that we can derive exactly the same equation 
by means of a skeleton equation for the order parame- 
ter. This means that the gaplessness of the BA mode is 
in fact a natural consequence of our truncation scheme. 
Together with the rescaled Ward identity (fl8|) and the 
gapless condition (fl9|) . the flow equations (fT5 |) -(fT7 |l form 
a closed system of RG equations for the five parameters 
(x), A, jl = Z{\i — £), r\ and 7 which can be solved 
numerically without further approximation. All bosonic 
polarizations can be expressed in terms of these parame- 
ters via skeleton approximations. For every value of the 
cutoff A this requires the numerical evaluation of a three- 
dimensional integral. The RG flow of the three-legged 
vertex 7 and the wave function renormalization factor 
Z at the unitary point in three dimensions is shown in 
Fig. H Moreover, in Fig. [3] we present our results for the 
order parameter (x), the single-particle gap A, and the 
chemical potential [i in units of the Fermi energy e F of 
a noninteracting system having exactly the same density 
as our flowing system. Note that cf is determined by 
the true density of the system which we calculate via the 
normal component of the one-particle Green function for 
a given value of the cutoff A. At the end of the flow we 
find for the renormalized quantities 

Ai/e F = 0.32, A/e F = G.61, (x}/cf = 0.59. (20) 

At the unitary point we may calculate the ground state 
energy per particle from Eq — 3/i/5 = 0.19. The above 
numbers should be compared with the mean-field results 
n/e F = 0.59 and A/e F = (x)/e F = 0.69. 



2 

A/e F)0 



FIG. 3: (Color online) RG flow the order parameter (x), the 
single-particle gap Aa, and the chemical potential fi at the 
unitary point in D = 3. Here, £f is the Fermi energy of 
a noninteracting system which has the same density as our 
interacting system. 



Our value [ijtp = 0.32 is smaller than the Monte Carlo 
results 0.44 (Ref. [IH) and 0.42 (Ref. 12l]), but it is quite 
close to the value 0.36 obtained by Haussmann et aZ. [l6| 
and agrees perfectly with the experiment by Bartenstein 
et al. 6j . Our value for the renormalized single-particle 
gap is within the error bars of the Monte Carlo simu- 
lations 11, [l2l | , a nd agrees with the FRG calculation 
by Diehl et aZ.jl3], while the FRG result for A/e F by 
Krippa [3] is only 3% smaller than the mean-field re- 
sult. In contrast to our approach, the FRG calculations 



of Refs. [3, El are based on a truncated gradient ex- 
pansion and do not distinguish between the quasi-particle 
gap and the order parameter, which are conceptually dif- 
ferent quantities. 

In summary, using a new truncation strategy of the 
FRG flow equations for the BCS-BEC crossover we have 
calculated the chemical potential, the single-particle gap, 
and the order parameter at the unitary point in three di- 
mensions and obtained reasonable agreement with exper- 
iments and other calculations. In contrast to the trun- 
cated derivative expansion of FRG flow equation for the 
BCS-BEC crossover used in Refs. 0, [ioj], our strategy 
is based on a truncation of the vertex expansion of the 
partially bosonized theory using Dyson-Schwinger equa- 
tions and Ward identities. Moreover, by introducing a 
momentum cutoff only in the bosonic sector, we can di- 
rectly calculate fluctuation corrections to the mean-field 
approximation, which serves as the initial condition for 
the FRG flow. The strong vertex correction of almost 
50% shown in Fig. [2] shows that Eliashberg type approx- 
imations are not quantitatively accurate close to the uni- 
tary point. 

This work was supported by the DFG via SFB/TRR49 
and by a DAAD/CAPES PROBRAL grant. 



5 



[1] D. M. Eagles, Phys. Rev. 186, 456 (1969). 

[2] A. J. Leggett, in Modern Trends in the Theory of Con- 
densed Matter, edited by A. Pekalski and R. Przystawa, 
Lecture Notes in Physics 115 (Springer, Berlin, 1980). 

[3] P. Nozieres and S. Schmitt-Rink, J. Low Temp. Phys. 59, 
195 (1985). 

[4] M. Drechsler and W. Zwerger, Ann. Phys. (Leipzig) 1, 
15 (1992). 

[5] M. Randeria, in Bose-Einstein Condensation, edited by 
A. Griffin, D. Snorke, and S. Stringari, (Cambridge Uni- 
versity Press, Cambridge, 1995). 

[6] M. Bartenstein, A. Altmeyer, S. Riedl, S. Jochim, C. 
Chin, J. H. Denschlag, and R. Grimm, Phys. Rev. Lett. 
92, 120401 (2004). 

[7] T. Bourdel, L. Khaykovich, J. Cubizolles, J. Zhang, F. 
Chevy, M. Teichmann, L. Tarruell, S. Kokkelmans, and 
C. Salomon, Phys. Rev. Lett. 93, 050401 (2004). 

[8] J. Kinast, A. Turlapov, J. E. Thomas, Q. Chen, J. Stajic, 
and K. Levin, Science 307, 1296 (2005). 

[9] G. B. Partridge, W. Li, R. I. Kamar, Y. Liao, and R. G. 
Hulet, Science 311, 503 (2006). 
[10] I. Bloch, J. Dalibard, and W. Zwerger, Rev. Mod. Phys. 
80, 885 (2008). 

[11] J. Carlson, S. Y. Chang, V. R. Pandharipande, and K. 

E. Schmidt, Phys. Rev. Lett. 91, 050401 (2003); S. Y. 

Chang, V. R. Pandharipande, J. Carlson, and K. E. 

Schmidt, Phys. Rev. A 70, 043602 (2004). 
[12] G. E. Astrakharchik, J. Boronat, J. Casulleras, and S. 

Giorgini, Phys. Rev. Lett. 93, 200404 (2004). 
[13] Y. Nishida, Y. and D. T. Son, Phys. Rev. Lett. 97, 

050403 (2006); Phys. Rev. A 75, 063617 (2007). 
[14] M. Y. Veillette, D. E. Sheehy, and L. Radzihovsky, Phys. 

Rev. A 75, 043614 (2007). 



[15] P. Nikolic and S. Sachdev, Phys. Rev. A 75, 033608 
(2007). 

[16] R. Haussmann, W. Rantner, S. Cerrito, and W. Zwerger, 

Phys. Rev. A 75, 023610 (2007). 
[17] R. B. Diener, R. Sensarma, and M. Randeria, Phys. Rev. 

A 77, 023626 (2008). 

[18] B. Krippa, larXiv:0704.3984V 3. 

[19] S. Diehl, H. Gies, J. M. Pawlowski, and C. Wetterich, 
Phys. Rev. A 76, 021602(R) (2007); S. Diehl, S. Flo- 
erchinger, H. Gies, J. M. Pawlowski, and C. Wetterich, 
larXiv:0907.2193l 

[20] S. Floerchinger, M. Scherer, S. Diehl, and C. Wetterich, 
Phys. Rev. B 78, 174528 (2008). 

[21] P. Strack, R. Gersch, and W. Metzner, Phys. Rev. B 78, 
014522 (2008). 

[22] We denote by K = (fc, iu>) and P = (p, iu) collective la- 
bels for momenta and Matsubara frequencies. Integration 
symbols are J K = ((3V)- 1 £ fc>w and f p = (/3V)' 1 £ p ^ 
where V is the volume and /3 is the inverse temperature. 
The corresponding normalization of the delta-symbols is 
&k,k' = l3Vo~u,u'&k,k' ■ We focus on the zero-temperature 
limit and infinite volume limit throughout this work. 

[23] C. Wetterich, Phys. Lett. B 301, 90 (1993). 

[24] T. R. Morris, Int. J. Mod. Phys. A 9, 2411 (1994). 

[25] F. Schiitz, L. Bartosch, and P. Kopietz, Phys. Rev. B 
72, 035107 (2005); F. Schiitz and P. Kopietz, J. Phys. 
A: Math. Gen. 39, 8205 (2006); P. Kopietz, L. Bartosch, 
and F. Schiitz, Lectures on the Renormalization Group - 
from the Foundations to the Functional Renormalization 
Group, (Springer, Berlin, to appear 2009). 

[26] See, for example, J. R. Schrieffer, Theory of Supercon- 
ductivity, (Benjamin, Reading, revised printing 1983). 

[27] L. Bartosch, H. Freire, J. J. Ramos Cardenas, and P. 
Kopietz, J. Phys.: Condens. Matter 21, 305602 (2009). 



